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FOREWORD 


The  purpose  of  the  program  is  to  investigate  cage  and  bearing 
instability  in  despun  antenna  bearings  due  to  changes  in  lubricant  properties 
The  program  is  being  conducted  at  Battelle's  Columbus  Laboratories,  505  King 
Avenue,  Columbus,  Ohio  43201,  for  the  Air  Force  Material  Laboratory,  under 
Contract  F-33615-74-C-5012.  (Project-Task  Number  7343  03  16.) 

Mr.  R.  J.  Benzing  and  Dr.  Mario  Rivera  (AFML/MBT)  of  the  Lubri- 
cants and  Tribology  Branch,  Air  Force  Materials  Laboratory  are  the  Project 
Engineers.  This  report  is  the  first  summary  report  under  the  contract  and 
covers  the  period  February  1,  1§74,  to  January  31,  1975. 

Mr.  J.  W.  Kannel  has  been  the  Principal  Investigator  and  has  had 
the  overall  responsibility  of  the  program.  He  has  been  assisted  at  Battelle 
by  Dr.  S.  S.  Bupara.  Most  of  the  experimental  work  published  herein  has 
been  conducted  by  COMSAT  Laboratories,  at  Clarksburg,  Maryland,  with 
Mr.  Chester  Pentlicki  acting  as  subcontractor  Project  Engineer. 

The  report  was  submitted  by  the  authors  January  31,  1975. 
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SUMMARY  AND  CONCLUSIONS 


i 


The  objective  of  the  program  has  been  to  perform  analytical  and 
experimental  investigations  of  cage  and  bearing  instability  in  low-speed , 
lightly-loaded  bearings  typified  by  bearing  for  despun  mechanical  syster,  s 
for  satellite  attitude  control.  The  analyses  are  being  conducted  at  Battelle' 
Columbus  Laboratories  (BCL).  They  involve  the  use  of  a computer  model  for 
the  bearing  which  predicts  cage  motions  as  a function  of  lubricant  properties 
and  bearing  design.  Empirical  inputs  of  lubricant  and  cage-spring  mass 
properties  are  used  in  the  analyses  to  ensure  accuracy  of  the  modeling-.  In 
addition,  experiments  are  being  conducted  at  COMSAT  Laboratories  to  observe 
cage  motion  as  a function  of  design  and  operating  conditions. 

The  analyses  have  involved  the  development  of  a specialized  version 
of  BCL's  bearing  dynamics  program,  BASDAP,  designated  as  BASDAP  II.  BASDAP  II 
is  distinguished  from  BASDAP  I in  that  it  uses  a series  of  closed  form 
solutions  for  in-plane  equations  of  motion  of  the  cage.  In  the  operation  of 
BASDAP  II,  the  motion  of  the  bearing  cage  is  traced  as  a function  of  time. 
Whenever  the  cage  impacts  a ball  or  race,  the  cage  momentum  is  altered  and 
the  cage  undergoes  a rebounding  action.  The  extent  of  rebounding,  and  hence 
a measure  of  the  energy  imparted  to  the  cage,  depends  on  the  extent  of 
elastohydrodynamic  lubrication  of  the  ball  race  interface,  the  viscosity 
of  the  bearing  lubricant,  and  the  spring  rate  properties  of  the  cage  material. 
An  analysis  of  the  rebounding  phenomena  has  led  to  the  development  of  a 
stability  criterion  (Equation  3). 

Two  types  of  experimental  evaluations  of  the  BASDAP  II  model  are 
being  pursued.  One  type  involves  a qualitative  correlation  between  actual 
spacecraft  observations  and  predictions  obtained  from  BASDAP  II.  The 
results  indicate,  both  from  the  stability  predictions  and  the  experimental 
observations,  that  typical  despun  bearing-lubrication  systems  operate  very 
near  the  threshhold  of  instability;  moreover,  relatively  small  changes  in 
the  external  conditions,  such  as  fluid  viscosity,  can  trigger  the  instability. 

The  second  type  of  experimental  validation  involves  operating 
bearings  under  various  conditions  of  lubricant  quantity,  lubricant  viscosity 
and  geometric  factors  and  of  monitoring  stability.  Experiments  conducted 
thus  far  with  Apiezon-C  at  room  temperature  and  a standard  flight  quality 
bearing  indicate  the  following. 

(1)  Stability  improves  with  increasing  lubricant 
quantity  for  near  starvation  condition  as  would 
be  predicted  by  BASDAP  II. 

(2)  Torque  variations  are  coupled  with  variations 
in  film  thickness  as  would  be  predicted  by 
BASDAP  II.  However,  it  is  not  certain  at  this 
stage  which  occurs  first. 

(3)  Some  instabilities  still  occur  for  adequate 
lubrication  conditions  as  could  be  expected  by 
virture  of  the  thin  EHD  film  and  high  viscosity 
of  the  lubricant. 
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INTRODUCTION 


Despun  mechanical  assemblies  (DMA's)  are  expected  to  function 
without  relubrication  or  maintenance  for  lifetimes  of  10  to  15  years  in 
space.  In  order  to  define  bearing  performance  and  reliability  in  such 
systems,  it  is  necessary  to  understand  the  failure  mechanisms  that  may  result 
in  loss  of  the  ability  of  the  system  to  perform  its  mission.  Furthermore, 
the  fundamental  failure  mechanisms  must  be  understood  in  order  that 
acceierated  tests  can  be  developed  to  evaluate  performance  life.  One  type 
of  failure  that  has  been  observed  in  space  for  DMA's  is  cage  (separator) 
instability  in  the  ball  bearings  which  has  caused  abnormal  torque  variation  u- 
and  serious  pointing  errors  in  the  antenna.  In  some  instances,  these 
instabilities  have  occurred  shortly  after  a change  in  the  thermal  environ- 
ment of  the  bearing  has  been  experienced.  Therefore,  serious  concern 
exists  about  the  useful  life  of  the  satellite  if  continued  deterioriation 
of  bearing  performance  should  occur.  In  addition,  it  is  possible  that 
increasingly  stringent  pointing  requirements  will  further  increase  the 
demands  on  the  DMA  bearing. 

In  order  to  develop  adequate  accelerated  life  tests  for  bearings, 

“ l\7,C*STy  t°  understand  aH  of  the  possible  failure  mechanisms  (cage 
nstability  has  been  shown  to  be  an  important  failure  mode)  so  that  separate 

Unfort.rnar  1tGS^  ^ 3nd  Perforraad  for  each  failure  mechanism. 

a . thj  State  of  the  art  is  such  that  an  understanding  of  how 

cage  instability  depends  on  bearing  design  factors  and  lubricant  properties 
is  not  available  as  yet.  In  the  absence  of  such  an  understanding,  it  is 
neither  possible  to  design  DMA  bearings  intentionally  to  be  stable,  nor  to 
devise  accurate  life  tests  to  evaluate  stability.  The  objective  of  this 
research  program  is  to  perform  analytical  and  experimental  investigations 
which  are  expected  to  enhance  this  understanding. 

The  analyses  are  being  conducted  at  Battelle’s  Columbus  Laboratories 
,hey  invo*ve  the  use  of  a computer  model  for  the  bearing  which  predicts 
cage  motion  as  a function  of  lubricant  properties  and  bearing  design;  this 
program  is  based  on  the  computer  program  BASDAP (1*2)  developed  earlier  hv 
BCL  and  is  called  BASDAP  II.  In  addition  to  the  analyses,  experiments  are 
being  conducted  at  COMSAT  Laboratories  to  actually  measure  cage  motion  as 
a function  of  design  and  operating  conditions. 


SIGNIFICANCE  OF  REPORTED  RESEARCH 


The  objective  of  the  research  program  reported  herein  has  been  and 
is  to  advance  the  understanding  of  cage  and  bearing  instability  in  despun 
mechanical  assemblies,  DMA's,  to  the  extent  that  life  test  methodology  can 
be  devised.  An  attractive  additional  objective  is  to  develop  design  criteria 
which  will  result  in  the  elimination  of  such  instabilities.  This  document 
represents  the  first  of  two  annual  reports  scheduled  for  the  contract,  and 
summarizes  the  groundwork  performed  toward  achieving  the  full  objectives 
of  the  program.  The  groundwork  has  been  oriented  in  such  a way  that  the 
remaining  research  efforts  will  represent  practical  application  of  the 
technology. 

In  reality,  significant  strides  have  been  made  toward  establishing 
not  only  a computer  program  for  studying  cage  instability  and  associated 
experimental  backup  but  also  toward  establishing  a stability  criterion  for 
the  bearing.  This  criterion,  if  validated  in  subsequent  work,  represents 
a momentus  step  forward  in  understanding  not  only  DMA  bearings  but  also 
most  types  of  angular  contact  bearings.  The  equation  for  the  stability 
criterion  can  be  written  (Equation  3 of  the  text). 
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where  is  the  effective  cage  mass,  Cg  is  the  ball-pocket  spring  constant,  and 


r = 2A  [eYPh(YP>,  - 1)  + 11 

^ h 1 ^Ph>2 


Here  A is  the  ball  race  contact  area,  h is  the  ball  race  film  thickness  and 
p^eYPh  is  the  viscosity  of  the  lubricant  between  ball  and  race. 

If  D is  greater  than  1,  the  cage  is  unstable.  Conversely,  if 
D is  less  than  1,  the  cage  is  stable.  The  bearing  dynamics  calculations, 
tne  laboratory  experiments,  and  spacecraft  observations  tend  to  substantiate 
this  criterion  for  a bearing  with  a ball-guided  cage.  However,  it  must  be 
noted  that  the  criterion  is  still  a preliminary  one. 

Assuming  the  validity  of  the  stability  criterion,  from  the  point 
of  view  of  bearing-lubricant  design,  the  DMA  designer  has  several  options. 

He  can  work  toward  a large  value  of  h,  a small  value  of  or  large  value 
of  Me  or  CSp  All  of  these  parameters  are  reasonable  design  variables. 

By  considering  the  criterion  from  the  point  of  view  of  life  test  methodology, 
the  test  engineer  can  accelerate  the  anticipated  lubricant  depletion  and 
attendant  starvation  in  the  bearing  to  force  premature  instability.  Likewise, 
he  can  evaluate  the  effect  of  cage  wear  on  the  spring  constant.  Also,  the 
role  of  wear  debris  on  viscosity  is  a possible  "failure"  related  parameter. 

In  general,  a stability  criterion,  if  it  holds  up  in  subsequent  research,  allows 
for  the  development  of  a model  law  for  bearing  failure. 


* M = 1/2  the  total  cage  mass. 
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PROJECT  ACTIVITIES 


The  bearings  used  in  DMA's  typically  are  very  light,  large  bore 
(‘100  rnra),  slow  speed  (-100  rpm),  lightly  loaded  (-100  lb)  bearings,  in  which 
normally  very  viscous  lubricants  are  used  to  ensure  adequate  ball— race 
lubricant  films.  The  cage  in  the  bearing  can  be  a conventional  race-guided 
configuration  or  a ball-guided  configuration  of  the  type  shown  in  Figure  1. 
Here,  a cone  surface  is  located,  inboard  of  the  normal  cage  cylinder,  to 
centralize  the  cage  location.  However,  since  the  cage  is  not  physically 
attached  to  any  other  bearing  element,  its  motions  are  subject  to  many 
uncertainties.  In  many  instances,  the  cage  will  incur  an  instability  (known 
as  a groan)  which  causes  serious  torque  fluctuations  in  the  bearing  apd 
pointing  inaccuracies  in  the  DMA.  The  purpose  of  the  research  presented 
here  has  been  to  elucidate  these  instabilities  so  that  test  methodology  can 
be  developed  to  evaluate  long-term  bearing  performance. 


The  specifications  used  for  the  antenna  bearing  study  are  given 
in  Table  I.  These  specifications  were  selected  by  mutual  agreement  between 
BCL,  AFML,  and  COMSAT.  This  bearing  represents  a typical  bearing  for  a DMA 
and  is  a bearing  with  a known  propensity  for  instability. 


DESCRIPTION  OF  BASDAP  II 


The  computer  program  BASDAP  represents  a very  complex  and 
generalized  program  for  the  evaluation  of  cage  dynamics  in  a bearing.  In 
the  BASDAP  computer  program,  the  cage,  itself,  is  treated  as  a 6 degree  of 
freedom  rigid  rotator.  This  computer  program  has  been  used  in  a variety  of 
research  programs  Including  antenna  bearing  studies.  However,  the  computing- 
time requirements  to  obtain  a convergent  solution  are  quite  long  and  some 
improvements  in  the  numerical  scheme  are  in  order  to  make  it  a more  useful 
tool  for  detailed  bearing  dynamics  studies. 


The  numerical  technique  used  in  BASDAP  involved  a fourth  order 
Runga-Kutta  numerical  scheme  for  the  integration  of  the  equations  of  motion 
for  the  cage.  This  involves,  then,  a point  by  point  numerical  integration 
of  the  equations  and  requires  time  steps  on  the  order  of  10“^  - 10-?  sec 
increments  to  obtain  convergence.  The  basic  approach  for  simplifying  the 
program,  is  to  obtain  closed  form  solutions  to  the  equations  for  at  least 
a portion  of  the  cage  motion.  This  simplified  computer  program  is  known 
as  BASDAP  II. 
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TABLE  I.  SPECIFICATIONS  USED  FOR  ANTENNA  BEARING  COMPUTER  STUDY 


Parameter 

Numerical  Value 

Number  of  balls 

16 

Ball  diameter,  in. 

0.59375 

Ball  mass,  lb 

0.022 

Separator  mass,  lb 

0.125 

Separator  pocket  diameter,  in. 

0.6130 

Outer-race  curvature,  percent 

52.5 

Inner-race  curvature,  percent 

52.0 

Bore  diameter,  in. 

3.543 

Pitch  diameter,  in. 

4.52765 

Contact  angle,  degrees 

26.14 

Cone  angle,  degrees 

67.5 

Axial  thrust  load,  lb 

60 

Shaft  speed,  rpm 

60 

f 


Approach  for  BASDAP  II 

Consider  the  bearing  of  Figure  1.  Assume  that  the  cage  starts 
in  motion  in  any  direction.  This  motion  will  continue  until  the  cage-pocket 
impacts  a ball  either  on  the  cone  or  on  the  pocket.  During  the  impact, 
the  rotational  motion  of  the  ball  will  impart  a momentum  to  the  cage  in  the 
direction  of  rotation.  In  addition,  under  some  conditions,  the  cage  will 
bounce  off  the  ball  so  that  the  resultant  motions  of  the  cage  will  tend  to 
be  away  from  and  radially  inward  (when  the  cage  hits  the  front  of  the  ball 
as  shown  in  Figure  1)  or  outwar I (when  the  cage  hits  the  back  of  the  ball). 

If  the  cage  were  race  guided  (i.e.,  if  it  contained  no  cone  surface),  then 
the  cage  could  also,  under  some  conditions,  impact  and  bounce  off  the  race. 

It  is  apparent,  then,  that  the  motions  of  the  cage  depend  heavily 
on  the  spring  characteristics  of  the  cage-ball  or  cage-race  interface,  on 
the  damping  of  the  motions  due  to  the  ball  slippage  on  the  races,  and  on  the 
friction  force  at  the  contact  interfaces.  A usable  analytical  model  for 
the  bearing,  then,  must  account  for  each  of  these  factors  in  a realistic 
manner.  And,  in  addition,  the  model  must  relate  these  factors  to  the  momentum 
transferred  to  the  cage.  With  a knowledge  of  these  momentum  transfer  forces, 
the  equations  of  motion  of  the  cage  can  be  solved. 

The  sequence  for  the  solution  to  the  cage  motions  is  as  follows. 

(1)  The  center  of  mass  motion  of  the  cage  (in  the  absence 

of  ball-cage  contact)  is  computed  by,  whenever  possible, 
a closed  form  solution. 

(2)  The  cage  rotational  motion  (in  the  absence  of  ball-cage 
contact)  is  computed  by  a closed  form  solution. 

(3)  The  ball-pocket  orientation  for  each  ball  is  scanned 
at  each  time  step  to  determine  the  onset  of  ball-to- 
pocket  contact. 

(4)  At  contact,  the  impact  equations  are  solved  by  a closed 
form  solution  and  the  impulse  transmitted  to  the  cage 

is  determined.  The  sequence  is  now  repeated  for  a large 
number  of  time  increments  and  the  forces  on  the  cage  and 
the  cage  dynamics  are  traced  as  a function  of  time. 
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Block  Diagram  of  BASDAP  II 


A block  diagram  for  the  BASDAP  II  computer  program  is  shown  in 
Figure  2 and  described  below.  Two  appraoches  are  used  for  Blocks  1 through 
3.  The  current  approach  is  to  compute  these  inputs  (exclusive  of  the  spring 
constants)  separately  and  to  input  them  into  BASDAP  II,  though  a more  general 
program  has  also  been  developed. 

Block  (1)  Input  Bearing  Design  Parameters 

This  block  sets  the  basic  parameters  for  the  bearing.  This  includes 
the  design  contact  angles,  the  bearing  speed,  the  lubricant  viscosity  including 
pressure-viscosity  effect,  and  the  bearing  geometry. 


Block  (2)  Statics  Analyses  of  Contacts 

(3) 

Using  the  theory  advanced  by  A.  B.  Jones  , the  contact  pressures 
between  balls  and  races  are  computed  as  are  the  contact  areas  and  contact 
angles.  These  analyses  are  described  in  Appendix  A and  represent  an  elasticity 
balance  and  allow  for  radial  loadings,  as  well  as  axial  loadings.  Centrifugal 
ball  loadings  are  also  inherent  in  the  computations  — though  they  are 
insignificant  for  a DMA  bearing.  Finally,  these  analyses  can  include  the 
computation  of  individual  ball  velocities  due  to  contact-angle  variations. 


Block  (3)  EHD  Films  and  Spring  Constants 

The  lubricant  film  at  the  ball-race  contacts  are  computed  using 
the  BCL  empirical  film-thickness  equation^).  These  calculations  utilize 
the  contact  stresses  and  velocities  developed  in  Block  (2).  Finally,  three 
constants  discussed  in  Appendix  B are  computed.  These  are 

C^  ~ The  ball-race  traction  constant 

Cg^  ~ The  ball-pocket  linearized  spring  constant 

Cg^  - The  cage-race  spring  constant  . 


Block  ( 4)  Initialize  the  Cage  Position  and  Velocities 

The  coordinate  system  for  the  cage  analysis  is  given  in  Figure  3. 
The  primary  coordinates  and  their  time-derivatives  for  the  in-plane  analyses 
are  as  follows 

p = The  displacement  of  the  center  of  mass  relative  to  the 
geometric  center  of  the  bearing 

g = The  angular  location  of  the  center  of  mass 

nt  = The  reference  angle  through  which  the  cage  has 
rotated . 
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FIGURE  2.  BLOCK  DIAGRAM  FOR  BASDAP  II  PROGRAM 
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is  set  initially  at  0. 


The  time  record  t^ 

Block  (5)  Time  Sequencing  for  Cage  Dynamics 

Time,  in  the  absence  of  ball-pocket  impact,  is  sequenced.  Whenever 
Block  (5)  is  entered  via  Blocks  (9)  or  (10),  the  incremental  time  sequence 
is  rezeroed.  However,  a continuous  record  of  the  total  time  is  maintained. 

Block  (6)  Computation  of  a,  a,  B,  B,  P,  and  p (No  Ball  Impact) 

CM  Radial  Motion  and  Displacement  (p  and  p).  The  center -of -mass 
of  the  cage  can  oscillate  radially  through  a maximum  of  2pm>  where  p is 
the  radial  clearance  between  the  cage  and  race,  without  contacting  the  race. 
For  the  case  where  p < p and  there  is  no  ball-cage  impact,  the.  only  external 
force  on  the  cage  is  that  due  to  centrifugal  loading  due  to  B . Conversely, 
for  a ball-guided  cage,  when  p > p , the  cage  contacts  the  race  and  a spring- 
type  force  is  applied  to  the  cage.  The  equation  of  motion  for  p,  then, 
depends  on  which  of  these  two  conditions  is  extant  at  any  particular  time. 

A special  subroutine  has  been  developed  which  computes  the  changes 
in  the  value  of  p and p through  any  time  increment.  The  analyses  for  this 
subroutine  are  given  in  Appendix  C. 

With  the  use  of  Equations  C-l  through  C-9,  it  is  possible  to  compute 
the  location  and  velocity  vectors  for  the  cage  CM  for  any  time  increment. 

In  the  operation  of  the  computer  program,  the  subroutine  is  supplied  values 
of  p.  and  p as  well  as  a series  of  functions  which  are  assumed  not  to 
varyiduring±the  time  increment  At.  The  computer  seeks  the  appropriate 
equation  (i.e.,  Equations  C-3  and  C-5)  and  the  time  that  this  solution  is 
active,  e.g.,  t . If  t is  less  than  At,  the  computer  switches  equations  at 
the  end  of  t and  repeats  the  logic.  This  progress  is  continued  until  the 
CM  is  traced  throughout  the  entire  At  increment.  The  method  for  tracing 
the  CM  motion  represents  a quasi  closed-form  technique  as  opposed  to  a 
finite-difference  solution,  since  it  utilizes  an  actual  solution  to  the 
motion  equation.  Thus,  the  subroutine  is  not  subject  to  the  same  numerical 
idiosyncrasies  inherent  in  finite-difference  techniques. 


Cage  Rotation  (a  and  a) (No  Impact).  In  the  absence  of  ball 
contact,  the  only  constraint  as  the  cage  rotates  is  associated  with  the 
friction  force  at  the  guiding  surface  when  p > pm.  In  the  computational 
techniques,  this  frictional  force  is  assumed  to  be  constant  throughout 
the  time  interval.  With  this  assumption,  it  can  be  shown  that 

la  = (Rc  - Rp)  Ff  , (1) 
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where  I is  the  mass  moment  of  inertia,  Y is  the  frictional  force  and'  (R  - R ) 
is  the  radial  location  of  the  f rictionalrforce  relative  to  the  center  ofCmassP 
(see  Figure  3).  Equation  1 can  be  readily  integrated  to  yield  a and  a as 
a function  of  time. 

B lock  (7)  Store  and  Print 

The  dynamics  computer  program  is  constructed  such  that  the  pertinent 
printouts  will  be  made  as  they  are  generated.  However,  to  minimize  the  amount 
of  output,  these  printouts  will  contain  only  the  pertinent  variables  such  as 
the  value  of  cage  location  vectors  (to  assess  stability)  and  the  ball-cage 
contact  forces.  In  addition,  the  data  are  stored  to  allow  for  direct  computer 
plotting. 


Block  (8)  Status  of  Ball-Pocket  Contact 

Each  ball  will  orbit  the  bearing  center  at  some  angular  velocity, 
“B(n)  (see  Figure  4),  dependent,  primarily,  on  the  contact  angles.  At  some 
condition,  the  angular  position  of  the  cage,  a,  will  be  such  that  ball- 
pocket  contact  will  occur.  The  condition  promoting  impact  is  described  in 
Appendix  D,  Equation  D-l.  This  impact  condition  is  related  to 

(a)  The  rotational  angle  of  the  cage 

(b)  The  center  of  mass  location  of  the  cage. 


Block  (9)  Momentum  Transfer  at  Ball-Cage  Impact 

As  discussed  in  the  approach  section,  a rebound  motion  of  the 
cage  will  occur  following  ball-cage  impact.  This  rebound  is  related  to 
the  spring-rate  of  the  ball-pocket  interface  and  the  viscous  resistance  of 
the  ball-race  interface.  In  simplified  form,  it  can  be  stated  that  the 
relation  between  the  rebound  velocity  and  the  impact  velocity  is  given  by 


where  u^  is  the  impact  velocity  of  a point  on  the  cage  with  the  ball  and  e 
is  a constant  of  restitution  which  is  related  to  spring  rate,  viscosity,  etc. 
Referring  to  Appendix  E,  the  value  of  e will  be  0 if  D <1,  otherwise  it 
will  be  the  value  given  by  Equation  E-6.  p 

In  order  to  determine  the  effect  of  impact  on  cage  motion,  it  is 
necessary  to  transform  the  cage  velocity  coordinates  (a,  0,  and  (3)  into 
effective  values  of  u.  and  v and,  following  impact,  it  is  necessary  to 

transform  the  value  or  u and  v back  into  cage  velocity  coordinates.  These 
c o o ° J 

transformations  are  given  in  Appendix  D,  and  consist  of  the  following  steps. 


n 
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(1)  The  value  of  8 and  p at  the  time  Impact  first  occurs 
are  transformed  into  rectangular  velocity  coordinate 
for  tbe  cage  CM,  x and  y (Equation  D-2). 

(2)  The  value  x and  y from  Step  1 along  with  the  initial 
value  of  a are  transformed  into  values  of  u and  v 
(Figure  4)  for  the  particular  pocket  incurring  impact 
(Equation  D-3). 

(3)  The  value  of  uq  is  now  computed  by  Equation  2 and  v 

is  computed  bv  Equation  D-9.  0 

(4)  The  values  uq  and  vq  are  transformed  into  the 
rectangular  velocity  components  x and  y for  the 
center  of  mass  (Equation  D-10) 

(5)  Finally,  the  values  of  x and  y are  transformed  into 
the  new  values  of  a,  8,  and  p. 


In  the  appendix,  the  prime  superscript  refers  to  the  CM  coordinates 
and  the  o and  i subscripts  refer  to  rebound  and  impact  conditions,  respectively. 
The  value  of  the  velocity  coordinates  d,  8,  and  p after  rebound  are  fed  into 
the  Block  5 of  the  computer  program  (Figure  2)  and  the  time  cycle  repeated 
until  a preset  number  of  total  time  step  have  been  incurred. 

A similar  procedure  is  used  if  ball-cone  contact  occurs.  The 
equations  for  this  type  of  contact  are  discussed  in  Appendix  F. 


COMPUTATIONS  USING  BASDAP  II 
Empirical  Inputs 


^he  success  of  any  analytical  modeling  of  a complex  system  such 

each  suiLbearin^/ePrdS  heavily  on  the  validity  of  the  specific  modeling  of 
each  subcomponent  in  the  system.  Two  inherent  assumptions  used  in  the  BASDAP 
program  are: 


(1)  The  ball  race  interface  is  modeled  as  a lubricated 
contact  where  the  tractions  are  modeled  by  bearing 
elastohydrodynamic  (EHD)  lubrication  theory. 

(2)  The  ball  cage  and  cage-race  interfaces  are  representative 
of  a dry  spring  interface  whence  the  normal  forces  are 
modeled  by  a spring  constant  and  the  traction  forces  are 
modeled  by  a coefficient  of  friction. 


In  order  then  to  ensure  realistic  BASDAP  calculations 
been  evaluated  experimentally  and  also  empirical  number 
for  inputs. 


these  models  have 
s have  been  used 
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Empirical  Evaluation  of  Ball-Race  Interface 


The  primary  factor  coatrolling  the  frictional  forces  at  the  ball 
race  contact  is  the  viscosity  of  the  lubricant  in  the  contact  zone  T^e 

csbaiC210  7 analysis  was  APiezon-C  (94.3  cs  at  100  F and  10.1 

, n / gaining  an  E.P.  additive.  A technique  using  a pair  of 
oiling  disks  was  used  for  the  rheology  evaluations  of  this  lubricant. 

in  Figure  5^  (°r  cylJnder)  viscometer  is  shown  schematically 

i j j T 1-inch  diameter  disks  (containing  a 16-inch  crown  radius! 
contact  to8ethf  lato  lubricated  (5000  rpm)  rolling  contact  to  the 
° 6lS  °f  interest  (100,000  to  200,000  psi).  Slippage 
introduced  between  the  disks  and  measurements  of  the  slip  rate  and  the 

use  ^ tractive  force  ia  »eaL“d  with 

use  of  a load  cell  which  constrains  the  upper  unit,  and  slippage  is 

etermined  by  measuring,  independently,  the  speed  of  each  disk  with 
electronic  counters  in  conjunction  with  magnetic  pickups.  With  the  apparatus 
traction  versus  slip  curves  are  obtained  for  various  WsTcSS^^e. 

Several  traction  curves  as  obtained  using  Apiezon-C  as  the  lubricant 

6;<  The  l0“-SliP  of  the  curved  ere 

sentially  linear,  which  implies  a linear  proportionality  between  the 
shear  stress  and  shear  rate  in  the  lubricant.  By  analyzing  this  Mnear 
portion,  inferences  about  lubricant  viscosity  can  be  made.  The  pressure- 
viscosity  curve  so  obtained  for  Apiezon-C  is  given  in  Figure  7 Under  the 
pressure  condition  of  interest  (>70,000  psi) , the  vi^Llly  of  this  ^brl 

laler  tMsyhighhvisf  'll  °f  10  CentiP°ise>  • Aa  will  be  discussed 

ter,  this  high  viscosity  can  cause  extreme  stiffness  at  the  ball-race 

junction  which,  in  turn,  can  cause  high  forces  to  occur  in  the  bearing 

v^c“lt5ed«a"4h£or''i1  Tt  Trf°a'  A1“'  ^ ln  FlgUre  7 are 

viscosity  data'  for  Vackote  obtained  by  standard  viscosity  measurement 


Empirical  Evaluation  of  Ball-Cage  Interface 

of  A"  anal>’tl“1  <P<pre8sio„,  based  on  theory 

evaluate^rhto’  I ball-cage  spring  rate,  Cs,  is  given  in  Appendix  B . To 
evaluate  this  spring  constant,  some  experiments  have  been  conducted  These 
experiments  involved  loading  a segment  of  the  cage  pocket  aglinst  a ball 
and  of  measuring  the  interface  deflection  characteristics.  8The  results  of 
these  experiments  are  given  in  Figure  8.  As  can  be  obse^ed  sulpllsing^ 

further ^olo  °d  £8  ^ '***  aChleVed  analytically . This  value  can  be 
further  improved  by  using  a value  of  Ec  of  0.8  x lo5  rather  than  105  psi. 

frictionalj?T^Sf-FLiCti0nal  FOrCe5-’  The  measurement  of  ball-cage 
segmellof  InTcl  i dlsCUSSed  in  Reference  2.  In  those  measurements,  a 
region  with  I Mn  h *'*V:atOT  was  loaded  ^inst  the  ball  in  the  pocket 
by  melnfof  foIr^J8  a?ra^*ment'  The  fiction  force  was  measured 

°ercen6d  7^  Phen°lla  Kparat^ha^ng  S«e“eJeirof  polosit^jo  7™ 

TheCIesultsPofCIhe’ball  ^ perCentJ’.as  suPPlied  by  COMSAT  Laboratories, 
results  of  the  ball-separator  friction  evaluations  are  summarized  in 
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FIGURE  8.  LOAD  DEFLECTION  CHARACTERISTICS  FOR  BALL-POCKET  INTERFACE 


Table  II.  For  the  cage  having  a 0.2  percent  porosity,  experiments 
were  conducted  w'th  and  without  oil  impregnation.  In  all  cases,  it  can  be 
observed  that  high  (and  nearly  constant)  coefficients  of  friction  were 
observed.  The  data  also  indicate  that  friction  in  the  presence  of  a 
lubricant  is  as  high  as  that  experienced  in  dry  contact.  The  values  ranged 
from  0.2  to  0.4  and  increased  with  increasing  porosity.  These  high  values 
for  friction  imply  that  the  attendant  large  normal  forces  between  ball  and 
separator  can  result  in  high  tangential  forces  which  might  affect  cage 
stability. 


TABLE  II.  DATA  FOR  BALL-CAGE  FRICTION 


(a)  Table  entry  is  the  dynamic  coefficient  of  friction. 


Preliminary  Computations  with  BASDAP  II 


Computer  Runs 


The  primary  efforts  thus  far  in  the  research  program  have  involved 
the  development  of  the  simplified  bearing  dynamics  program,  BASDAP  II  and 
preparation  of  empirical  inputs.  In  essence,  this  program  computes  in-plane 
cage  motion  using  basically  a closed  form  integration  scheme  for  the  equation 
of  motions.  The  program  can  be  used  to  evaluate  both  a ball  guided  and  a 
race  guided  cage.  In  addition,  two  specific  cases  were  evaluated  for  a 
ball  guided  configuration,  one  case  involved  using  a value  of  lubricant 
viscosity  that  produced  high  ball-race  tractions,  whereas  the  second  case 
involved  low  ball-race  tractions. 


Graphs  of  the  prediction  of  cage  oscillation,  a,  versus  time  are 
given  in  Figure  9 for  a short  period  of  bearing  operation.  At  0.0149  sec 
the  high  traction  case  produced  intolerable  a motions  which  aborted  the 
computations,  (i.e. , the  equations  predicted  cage  instability).  Conversely, 
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for  the  low  traction  case,  the  a diminished  toward  zero  after  about 
0.016  sec  of  real  time  indicating  definite  stability  of  the  predicted 
cage  notion.  Mechanistically,  then,  the  viscosity  of  the  lubricant  at 
the  ball-race  contact  is  essentially  responsible  for  the  damping  of  the 
cage  motions.  That  is,  high  viscosity  (i.e.,  high  tractions)  causes  the 
cage  to  bounce  off  the  ball  whereas  low  viscosity  causes  the  ball  to 
absorb  the  cage  motions. 

The  prediction  given  in  Figure  9 is  extremely  interesting.  If 
in  further  evaluations  they  are  found  to  be  realistic,  then  a real 
understanding  of  cage  instability  is  evolving.  As  discussed  in  the 
subsequent  section,  this  stability  criterion  should  be  subject  to  "model" 
laws  which  can  be  used  to  evaluate  cage  design  concepts  as  well  as  life 
test  methodology  for  cage  instability. 


Stability  Criterion 


One  advantage  of  the  type  of  computer  technique  used  for  BASDAP  II 
(as  opposed  to  the  Runga-Kutta  integration  of  the  equation  of  motion  used 
in  the  original  BASDAP),  is  that  it  is  possible  to  determine  what  aspect 
of  the  computer  model  is  triggering  an  instability.  The  two  cases  discussed 
in  the  preceding  section  involved  evaluations  of  low  and  high  ball-race 
tractions.  These  tractions  only  affect  the  BASDAP  II  calculations  through 
Equation  E-4  of  Appendix  E.  This  equation  appears  in  the  form 


D 

P 


(3) 


where  M is  effective  cage  mass,  Cg  is  the  ball-cage  spring  constant  (from 
Appendix  B and  the  preceding  sectioi)  and  C is  a viscous  term  for  the  ball 
race  interface  (also  discussed  in  Appendix  B). 


If  D is  greater  than  1,  the  cage  will  "bounce"  off  the  ball  and 
the  cage  will  Rave  a propensity  to  instability.  Conversely,  if  Dp  is  less 
than  1,  the  cage  motion  will  be  absorbed  by  the  ball  and  the  cage  will  tend 
to  be  stable. 


EMPIRICAL  EVALUATION  OF  BASDAP  II 


The  credibility  of  the  bearing  dynamics  model  depends  on 

(1)  How  well  the  computer  output  results  correlate 
with  actual  space-bearing  operation. 

(2)  How  well  the  output  data  correlates  with  the 
results  of  laboratory-controlled  experiments. 
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BASDAP  Evaluation  of  Spacecraft  Bearings 


Cage  instability  in  actual  DMA's  has  been  detected  by  accelero- 
meters attached  to  the  spacecraft.  During  such  instabilities  the  normal 
60-80  Hertz  wobble  noise  of  the  cage  deteriorates  into  an  audible  groan 
from  the  bearing.  Such  cage  behavior  has  been  seen  in  laboratory  tests, 
in  preflight  tests  of  spacecraft  hardware  and,  as  mentioned  above,  in  actual 
space  operation,  A cure  for  the  bearing  groaning,  in  space,  has  been  to 
increase  the  bearing  cavity  temperature  from  the  nominal  70  F to  around 
100  F.  One  check  of  the  bearing  dynamic  calculations  then  would  be  to 
compare  stability  prediction  for  the  two  temperature  situations.  This  was 
done  for  the  case  where  the  lubricant  is  Apiezon-C  and  the  bearing  is  the 
one  described  in  Table  I. 

Typical  condition  for  a spacecraft  bearing  for  the  two  temperatures 
are  summarized  in  Table  III.  As  evidenced  by  the  value  of  D , the  bearing 
conditions  are  very  close  to  the  stability  threshhold.  NotePthatD  can 
realistically  vary  widely  on  either  side  of  unity  depending  on  bearing 
conditions.  The  stability  criterion  for  BASDAP  II,  thus  appears  to  predict 
very  reasonable  conditions  for  a flight  bearing. 


TABLE  III.  STABILITY  ANALYSIS  FOR  FLIGHT  BEARINGS 


Parameter 

Units 

Case  A 

Case  B 

Temperature,  T 

F 

70 

120 

Lubricant  Base  Viscosity,  p 

cs 

250 

60 

Lubricant  Intercept  Viscosity, 
^i 

cp 

3.6  x 104 

0.9  x 104 

Pressure  Viscosity  Exp,  y. 

.-1 

psi 

3.15  0.10'5 

3.15  0.10'5 

Hertz  Pressure,  P, 
n 

psi 

78,000 

78,000 

Contact  Area,  A 

in-4  . 2 

10  in. 

1.73 

1.73 

Effective  Cage  Mass,  M 

e 

O 

lb  ’ sec  / in. 

3.2  x 10"4 

3.2  x 10'4 

EHD  Film  Thickness,  h 

10  ^ in. 

10.0 

4.0 

Ball-Cage  Spring  Rate,  C (linear] 

lb/in. 

2 x 103 

2 x 103 

1 

Ball-Race  Damping,  C 

0.54 

0.33 

Damper  Parameters,  Dp 

1.82  (Unstable) 

0.7  (Stable) 
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Laboratory  Evaluations  of  BASDAP 


Most  of  the  experiments  are  to  be  conducted  by  COMSAT  Laboratories  i 

using  flight  equivalent  bearings.  These  experiments  are  to  involve  measure-  ' 

ments  to  assess  the  stability  characteristics  of  the  bearing  under  various  i 

conditions  of  lubrication  and  for  various  type  of  configurations.  The 
measurements  will  include 

(1)  Dynamic  torque  fluctuations. 

(2)  Extent  of  V r-"e  EHD  lubrication  films  using  an 
electrical  conductivity  technique.  That  is,  a voltage 
is  imposed  across  the  beai ing  and  the  interruptions 

to  that  voltage  to  metallic  contact  through  the  EHD 
film  is  observed. 

(3)  Some  actual  cage  motions.  f 


Effect  of  Lubricant  Quantity  on  Stability 

Two  types  of  experimental  series  are  being  conducted  to  validate 
the  analysis.  One  experimental  series  is  being  conducted  to  evaluate  the 
propensity  of  the  cage  to  become  unstable  as  a function  of  lubricant 
quantity  in  the  bearing.  The  second  series  of  experiments  will  be  conducted 
using  a homologous  series  of  lubricants  (such  as  the  Kendall  SRG  series)  to 
establish  the  effect  of  viscosity  on  instability.  A minimum  of  two  lubri- 
can  will  be  used,  including  Vackote.  The  exact  experiments  will  depend 
on  tne  data  developed. 


Some  preliminary  data  have  been  generated  in  the  oil  quantity 
experiment.  The  experimental  sequence  for  these  tests  is  as  follows. 

(1)  The  bearing  is  operated  dry  and  the  torque  behavior 
measured  for  100  rpm  motion,  60-lb  load,  and  at  room 
temperature . 

(2)  Two  drops  of  lubricant  (Apiezon-C)  are  injected  into  the 
bearing  on  two  balls  180  degrees  apart.  Each  drop  weighs 
approximately  4 mg  and  is  a mixture  of  20  percent  oil 
and  80  percent  Freon.  Torque  and  EHD  films  are  then 
monitored  as  function  of  time  for  several  minutes. 

(3)  Additional  quantities  of  oil  are  added  to  the  bearing  and 
the  EHD  film  and  torque  behavior  noted. 


Typical  strip  chart  recordings  of  bearing  behavior  are  given  in 
figures  10  through  13.  (Charts  are  shown  in  increasing  degree  of  lubrica- 
tion as  outlined  in  1 through  3 above.)  It  can  be  observed  that  there  is 
a definite  tendency  of  the  bearings  to  stabilize  as  lubricant  is  added. 
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However,  there  is  still  some  evidence  of  instability  at  full  film  lubrication. 
This  is  consistent  with  the  results  shown  in  Table  III.  There  is  also  an 
experimental  relation  between  loss  of  EHD  film  and  torque  fluctuation  It 
is  not  clear,  however,  whether  the  loss  of  film  comes  first  or  the  instability. 

tability  has  also  been  observed  to  be  speed  dependent  occurring  most 
often  above  40  rpm. 


The  data  from  the  laboratory  experiment  conducted  thus  far  are  in 
general  agreement  with  the  stability  discussion  of  the  preceding  section 
and  with  Equation  3.  Specified  points  of  comparison  are  as  follows. 

(1) 

Stability  is  improved  with  lubrication  and  with 
lubricant  quantity  (i.e.,  h in  Equation  3). 

(2) 

Stability  is  directly  affected  by  film  thickness 
fluctuation  or  vise  versa. 

(3) 

Even  with  adequate  lubrication,  some  instabilities 
still  occur  which  is  consistent  with  the  calculation 
of  Table  III. 

There  are  still  some  areas  that  merit  further  analysis  such  as 
the  exact  cause  relationship  between  film  thickness  and  torque  fluctua- 
tion and  the  reason  for  the  critical  speed  for  instability.  These,  as  well 
as  other  factors,  will  be  explored  in  much  more  detail  both  experimentally 
as  well  as  analytically  in  the  second  half  of  the  research  program. 


Effect  of  Bearing  Design  Factors  on  Stability 

Several  experiments  are  being  planned  to  evaluate  the  role  of 
bearing  design  factors  on  cage  instability.  These  experiments  which  have 
just  been  initiated  will  include  evaluation  of  the  following  design  factors. 

The  desi8n  factors  being  considered  for  experimental  evaluations 


(1)  Contact  Angle.  Here  it  is  planned  that  one  contact 
angle  (in  addition  to  the  26  degree  angle)  will  be 
evaluated.  This  angle  will  be  significantly  different 
than  26  degrees,  although  the  exact  angle  is  dependent 
on  ball  procurement  by  COMSAT. 

(2)  Effect  of  Cage  Loading.  Experiments  will  be  conducted 
for  various  orientations  of  the  g-vector  relative  to 
the  spin  axis. 

(3)  Cage  Pocket  Configuration.  Experiments  will  be  conducted 
for  a ball-guided  cage  with  at  least  two  cone  angles. 

In  addition,  two  cage  pocket-materials  (phenolic  and 
metal)  will  be  used  to  evaluate  the  effect  of  cage 
stability  of  the  pocket  spring  rate. 
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(4)  Effect  cf  Cafe  Control  Surfaces.  Experiments  will  be 
conducted  with  a race-guided  cage  as  well  as  a ball 
guided  cage. 

ah  f pvnpr iments  will  be  conducted  at  a speed  of  100  rpm 

and  a load  of  60  lb  at  ambient  tempetatute.  The  lubricant  for ^^eteminLg 

The  above  bearing  test  program  is  quite  general  in  scope  and  is 

it  is  very  difficult  to  be  much  more  specific  in  planning  tne  e v 
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NOMENCLATURE  FOR  APPENDICES 


Half  length  of  ball  race  contact  area. 

Half  width  of  ball  race  contact  area. 

Subscript  refer  to  contact  at  cage  cone  (Appendix  F). 

Race  cage  spring  constant. 

Ball  cage  interface  spring  constant  (Equation  B-10). 

Ball  race  viscosity  damping  constant  (Equation  B-4) . 

Center  of  mass. 

Damping  constant  for  ball  race  interface  (Equation  E-l). 
Ball  diameter. 

Constant  of  restitution  of  ball  cage  impact.  (Equation  E-6) 
Friction  coefficient. 

Race  curvature. 

Refers  to  force. 

Ball-race  tractive  force. 

Centrifugal  force. 

EHD  film  thickness. 

Subscript  for  inner  race  (or  input  condition). 

(Equation  D-l) . 

Ball-race  spring  constants. 

Cage  mass. 

Effective  cage  mass  (1/2  M ) 

c ' 

Inner  to  outer  race  misalignment. 

Refers  to  nth  ball. 

Number  of  balls  in  bearing. 

Subscript  for  outer  race  (or  output  condition). 

Ball  loading. 

Maximum  Hertz  pressure. 

Ball  radius. 

Radius  of  Pitch. 

Any  normal  position  vector  or 
Subscript  refers  to  static  contact. 

T ime . 

Velocity  in  y direction. 

Velocity  in  y direction. 


NOMENCLATURE  FOR  APPENDICES  (Continued) 


v - Velocity  in  y direction, 
x - Coordinate  variable, 
y - Coordinate  variable. 

ag(n)  = Ball  location  angle  for  ball  n. 

ac (n)  ~ Cage  location  angle  at  ball  n. 

0 ~ Cage  center  of  mass  rotational  angle. 

Y - Pressure  coefficient  of  viscosity. 

^ - Deflection  of  an  interface. 

A»  Ao’  Ai  _ Bal1  race  interface  total  deflection 

e - Radial  clearance  between  ball  and  pocket. 

<p  - Cage  cone  angle  (Figure  F-l). 

0c,es’  0o,0i  " Bearing  contact  angle  (Appendix  A). 

p - Radial  location  of  cage  center  of  mass  (Figure  3). 
Pm  " Maximum  value  of  p . 
v - Poisson's  ratio. 


APPENDIX  A 


BEARING  STRESS  COMPUTATIONS 


n,  , „ T^e  PurP°se  of  this  Appendix  is  to  outline  the  approach  for  deter- 
mitig  the  ball-race  contact  stresses  and  dimensions.  This  approach 
involves  first  computing  the  stresses  for  a static  (nonrotating)  bearing 

and  fnpn  nnmr»i irinn  ~ r . r-  - ~ &/ 


• « ~ vHuiuuiaiiiie  j Dear  i 

and  then  computing  the  effect  of  centrifugal  forces  on  contact  angles. 
The  equations  have  been  abstracted  from  the  works  of  A.  B.  Jones. (3>7) 


Static  Analyses 


,,,  ,,  Dn^er  static  contact,  the  bearing  can  be  subject  to  both  axial 
(induding  misalignment)  and  radial  loads.  For  the  purpose  of  the  analysis 
this  implies  that  the  deflection  of  each  ball-race  interface  in  the  bearing 
must  produce  forces  which  sum  to  the  applied  radial  and  axial  loads  § 


Deflections 


Referring  to  Figmre  A-l,  the  vertical,  A , and  horizontal,  A,, 
separations  of  the  race  curvature  centers  can  be  expressed  as  1 


Aj(n)  = Bd  sin6  + 6 + m R cos  ( — (n-1)  V) 

u u P V N J 


(A-l) 


A2(n)  = Bd  cos  0rf  + 6^  cos  \ 

NB  ^ 


(A -2) 


Here  d is  the  ball  diameter,  0d  is  the  design  contact  angle,  m is  the 
salrgnment,  R is  the  pitch  radius,  n refers  to  the  nth  ball,  NR  is  the 

brdaet^dlLballS*  ^ 3nd  61  -e  the  radial  deflections  to 


B=  fo+fi  "I. 


where  f , are  the  outer  and  inner  race  curvatures.  (Bd  is  the  distance 
between  durvature  centers.) 


The  total  deflection  of  the  ball  race  interface  (for  each  bc’l) 
can  be  written  ' 


Mn)  = ^ Ax2  + A 22  - Bd 


and  the  contact  angle  after  deflection  becomes 


0 - t.„-1  ^ 

5 A2 


Wll 


Outer  Race 


Inner  Race 
Curvature 
Center 


Outer  Rac* / 
Curvature  C * 
k Center/' 


Design 


Inner  Race 


FIGURE  A-l . NOMENCLATURE  FOR  STATIC  ANALYSES 


Force  Equilibrium 

deflection^  at  bal1  contact  can  be  written  in  terms  of 

3 

ps(n)  = K(n)  A2  , 

where  K is  the  interfacial  spring  constants  given  by 


l 

K 


o i 


(A-7) 


race  contact.  Tfe^constintt^rl^def i^ed^for"8^?  f°r  the  °Uter  and  inner 
B-7,  and  B-8.  defined,  for  example,  in  Equations  B-6, 


Ignoring  gyroscopic  effects,  the  force 

NB 

fa  " E ps  sin  6S 

n=l  b 


summation  requires  that 
* (A-8) 


and 


nb 

E P 
n=l 


. cos  eo  sin 
s s 


retain  1 

L nd  J 


(A-9) 


to  A-9  have^been^olved  ^sinj  a^oible"3  “m*1  force  components.  Equations  A- 

(Erua^^^^A.r^rLfr^i-irt1;^  rV  r v?1 

guesses  of  «0  and  §1  are  made  until  FA’.nS  V™  t^pp^d  ^g 

Effect  of  Centrifugal  Force 

the  inner  2 ZLsT  ^ th“  force  niters 

affect  the  ball-to-ball  load  sharing!  ( ass°oiated  forces)  but  does  not 

Deflection  Equations 

contact  angle  s^in^of  "course*  ‘Zb  “"d  T'  1*“  “ «“  *■  the 

Figure  A-j'it  can  li  °JL?5£’  ^ S<’“”h‘'t  dlft“o»t.  Referring  to 


Ao(n) 


J xL2  + 


and 


X2  (f0  ‘ °-5>  d , (A-10)  • 

**(n)  = /^X,2  + (A2  - X22  - (fi  - 0.5)  d . (A-ll) 
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FIGURE  A-2.  NOMENCLATURE  FOR  CENTRIFUGAL  FORCE  EFFECTS 


Likewise 


3 3 

P1  - KtAt  and  PQ  - KQ  A2  . 

(A-12) 

Also 

-1  X1 

60  ' tan  - , 

(A -13) 

and 

- 1 !“  ^ 1 " 

(A -14) 

Force  Equilibrium 

The  force  equilibrium  between  inner  and  outer  race 
gyroscopic  force  effects,  can  be  written 

ignoring 

p0  Sin  eo  = pi  sin  6i 

(A-15) 

and 

P0  cos  0O  = pi  cos  6i  + Fc 

(A -16) 

The  assumption  that  centrifugal  effects  do 
sharing  implies  that 


not  affect  ball-ball  load 


P.  sin  0.  = P sin  0 
i i s s 


(A-17) 
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Solution 


By  combining  Equations  A-15  and  A-16  there  results 

F 

_1 = — 1 £ (A-18) 

tan  6 tan  0,  Po  sin  6 

o i s s 


By  guessing  (by  a nesting  process),  a value  of  0^,  9Q  can  be  computed.  With 
these  values  of  0^  and  0O,  and  X2  can  be  computed  from  Equations  A-13 
and  A-14  for  the  known  (from  the  static  analysis)  value  of  A^  and  A2. 
Equations  A-10  to  A-12  can  now  be  used  to  determine  a value  of  Pj>.  The 
correct  guess  of  0^  will  produce  a value  of  P^  which  is  consistent  with 
Equation  A-17. 


I 
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APPENDIX  B 


1 -J 
( 


SPRING  AND  DAMPING  CONSTANTS 


Ball-Race  Force -Traction  Relationship 


If  it  is  assumed  that  the  ball  and  race  in  a bearing  are  separated  by 
a thin  layer  of  lubricant  of  thickness,  h,  then  the  ball-race  tractive  force 
can  be  written 


II  T dA 

A 


(B-l) 


where  T is  the  shear  stress  and  A is  the  contact  area.  It  can  be  shownv  J that 


(5) 


T - - HeVP  f 


By 


(B-2) 


Here  AU  is  the  slip  at  the  ball-race  interface,  is  the  base  viscosity,  p is 
the  contact  pressure,  and  y is  a constant  pressure  viscosity  coefficient.  It 
is  assumed  that  the  contact  pressure  distribution  is  Hertzian,  i.e., 


2 2 
x z_ 

V ■ Ph  V1  ' ~2  ■ 2 

a b 


(B-3  ) 


where  p^  is  the  maximum  contact  pressure  and  a and  b are  the  major  and  minor 
axis  of  the  contact  ellipse. 


By  combining  Equations  B-l  through  B-3  there  results 


(5) 


F = C AU 
T p 


2TTabpi  vp. 

°u " rr~2  >e 

hy  p. 


1)  + 1] 


(B-4) 


For  a given  ball  race,  contact  situation  the  ball  race  traction  is 
then  dependent  only  on  a constant,  C^,  and  the  slip  condition. 
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.Normal  Force  Between  Ball  and  Cage  Pocket: 


A reasonable  assumption  for  the  ball-cage  contact  is  that  it  represents 
an  elastic-spring  mass  system.  The  force  deflection  relationship  for  two  elastic 
bodies  in  contact  can  be  written  (from  Sealy  and  Smith(8),  page  364). 


r-  "C.  2 2 i 

LC6(A+B)J  a " s 6 * (B-5) 

where  P is  load,  and  C are  constants,  and 


A = 


(B-6) 


A 


1 

A + B 


1-v. 


(B-7) 


Ilere  6 is  the  deflection,  r is  radius,  v is  Poisson's  ratio,  E is  Young's  modulus 
and  the  subscripts  p and  B refer  to  pocket  and  ball,  respectively. 

Similar  equations  can  be  written  for  the  cage-race  contact.  Here  rs 
would  be  the  race  radius  and  r would  be  the  cage  radius.  The  parameter  B would 
also  be  somewhat  modified  to  accommodate  the  different  contact  geometry  condition. 

The  constant  and  Cg  can  be  approximated  by 


Cb~ 


•635  (t 


-0.1722 


and  C ss  3 


214  I? 


■0.471 


(B-8) 


Equation  B-5  involves  the  use  of  a single  constant,  Cs,  and  the  deflection  to 
determine  ball-cage  forces. 


A very  useful  simplification  would  occur  if  Equation  B-5  could  be 
approximated  by  a linear,  rather  than  a 3/2  power  spring  in  the  form 


P C 6 
S1 


(B  —9  ) 


If  a typical  value  of  P is  say  0.2  lb,  then  Equations  B-5  and  B-9  would  produce 
the  same  load  for  this  condition  provided  that 

2 

Cs1W°*6Cs3  * (B-10) 

Similarly,  can  be  computed  for  a cage-race  contact. 
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APPENDIX  C 


CENTER  OF  MASS  LOCATION  (INCLUDING  RACE  GUIDED  CAGE  INFLUENCE) 


Center-of -Mass  Location 


Referring  to  Figure  3,  the  equation  of  motion  for  radial  displacement 
of  the  cage  CM  can  be  written  (in  the  absence  of  ball  contact) 

P + M1  " Fc  = 0 ’ (C_1) 

c 

where  M^  is  the  mass  of  the  cage,  p is  the  radial  displacement  and 


F = 
s 


and,  it  is  helpful  to  define 


F = 
c 


Here  pm  is  the  maximum  free  movement  of  the  cage  and  Cg  is  an  assumed  linear 
spring  constant.  Thus,  Fg  characterizes  the  cage-race  radial  contact  force 
and  Fc  characterizes  the  centrifugal  force  acting  on  the  cage  CM. 


Solution  for  p < pm 


For  this  condition  the  solution  to  Equation  C-l  can  be  written 


C.  sinh  (A  + Pt)  if  C.  real 
11  . (C-3 ) 

iC^  cosh  (A  + Pt)  if  C^  imaginery 


r Cs  (p  - Pm)  if  p > p 

l 0 


m'  “ r ^m 
if  p < p 

m 


(C-2) 


Fc,  such  that 


e2p 

if  p < pm 

P2pm 

If  p > p_ 
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By  designating  and  as  the  initial  radial  location  and  velocity  of  the 
CM,  respectively,  then 


and 


■{ 


-1  , 

sinh  1 

V.  cx  y 

if  Cx  real 

-1  , 

r pi  n 

• 

cosh  1 

v.  ±c1  y 

1 if  imaginary 

(C-4) 


Equation  C-3  can  be  inverted  to  determine  the  time  when,  for  example,  p *■=  pm 
as  well  as  being  used  to  compute  p for  a given  time  step. 

Solution  for  p > pm 


The  solution  of  this  equation  can  be  written 


p = C'  sin  (A  ' + Bt)  + C2  , 


(C-5  ) 


where 


if  and 

c 


C2  " C1  + ^ ^ ’ 


(C-6) 


Using  the  same  initial  conditions  as  above,  it  can  be  shown  that 


_ i r ® ^ p a 

A ' = tan  f — 

L P_. 


B(P,  " C2)  ^ 


and 


(C-7 ) 


:1  = BcosTa7) 


(C-8) 


APPENDIX  D 


MOMENTUM  TRANSFER  TO  CAGE  CYLINDER  AT  IMPACT 


Momentum  Transfer  at  Ball— Cage  Impact 


the 


For  each  time  step,  the  computer  program  determines  value  of 
ollowing  parameters  and  their  time  derivative  (see  Figure  D-l) 


°c (n)  = The  angular  position  of  the  cage  at  any  ball  location-n 

Vn)  * The  angular  position  of  the  ball-n 

0 = The  angular  position  of  the  cage-CM 

p = The  radial  position  of  the  cage-CM. 


Impact  Criterion 


Impact  will  occur  between  ball  and  cage  for 
less  than  zero,  where 


any  ball-n  when  CHK 


is 


CHK 


r [»c  + p/Rp  sin  (P  - Q'c ) ] - c*B  + e 
“b  “ fac  + p/RP  sin  “ «c) ] + e 


front  (ic  * -1) 
back  (ic  * 1). 


(D-l) 


RP  Is,  JhE  Pltch  radl“s  and  £ ^ the  radial  clearance  between  ball  and 
pocket.  At  front  contact  (see  Figure  D-l)  i is  set  equal  to  (-1)  and  at 
back  contact  i ^ is  set  equal  to  +1.  c H 1 ' aC 


Impact  Velocity 

At  impact,  the  components  of  velocity  of  the  CM  of  the 
written  (in  terms  of  known  parameters)  as  follows 


cage  can  be 


X£  = sinP  + pp^  cosP 
y-  = Pj  cosP  - pg  sinP 

where  the  i subscript  refers  to  impact  condition. 


The  approach  velocity  of  the  cage  to  ball,  u,  can  be  written 

Ui  = RP  <°c  ' V + *t  cosac  " y[  sinac>  and 

v * x'  sino  + y ' coso 
1 1 c i c 
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Momentum  Transfer 


During  impact,  the  velocity  component  u and  v are  altered.  The 
alteration  of  u is  as  described  in  Appendix  E.  Following  impact, 

Uo  “ “ui  e » (D-4) 

where  e is  a constant  of  restitution.  The  o and  i refer  to  output 
velocity  and  input  velocity,  respectively.  The  change  in  the  velcotity 
v can  be  derived  from  momentum  considerations  as  follow 


Mdv  = 

F dt 

y 

, since 

(D-5) 

Mdu  = 

F dt 

X 

, and  since 

(D-6) 

F = 

y 

fF 

X 

» 

(D-7) 

there  results 


Vo  ~ vi  = f (u0  “ . (D-8) 

where  f is  the  friction  coefficient  and  F and  F refer  to  force  components 
at  the  cage  impact.  In  reality,  however  v cannXt  exceed  the  linear  velocity 
of  the  ball  surface  vfi  so  that 

f [vi  + f(uo  ' ui)ic]  if  Vo  < VB  (D-9) 

v *=  j 

v v„  i otherwise. 

B c 


Dispart  Velocity 

The  alteration  of  the  velocity  components  u and  v will,  of  course, 
modify  the  center  of  mass  velocity  components  or* 


Cu  - u . v 

1 J cos orc  + (vq  - Vi)  sina.,  and 
K ■ y't = ’(  °2  ')  Sinac  + (vo  " V cosoc 


(D-10) 


Finally,  by  inverting  Equation  D-2,  there  results 


p = x1  sinP  + y'  cosP,  and 
o o 


• • • f 

pP  = x ' cosP  - y sing 


(D-ll) 


* 


2 

It  is  assumed  here  that  I = M Rp  t where  I is  the  cage  mass  moment  of 
inertia  and  is  the  cage  mass. 
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APPENDIX  E 


IMPACT  OF  BALL  AND  CAGE 


APPENDIX  E 


IMPACT  OF  BALL  AND  CAGE 


Under  impact  conditions,  a force  will  be  exerted  on  the  cage  which, 
if  sufficiently  large,  can  affect  the  cage  motion  adversely  and,  possibly, 
can  perpetrate  an  instability.  To  examine  the  condition  under  which  an 
instability  might  occur,  consider  the  force  diagram  of  Figure  E-lb.  Here 
cp  and  Cs  are  the  constants  described  by  Equations  B-4  and  B-9,  Me  is  the 
effective  mass  of  the  cage,  6 is  the  deflection  of  the  cage  at  the  ball- 
cage  interface,  and  s is  the  cage  velocity.  For  equilibrium  of  the  forces 
ignoring  ball  mass  in  comparison  to  ball  race  friction,  it  can  be  shown 
that 


C (6-s)  = - C 6 
\i  sx 

(E-l) 

-C  6 = M s 
s e 

where  M is  the  effective  mass  of  the  cage.  (This  can  be  shown  to  equal 
/ MC.J  By  performing  an  integration,  the  above  two  equations  can  be 
combined  to  yield  (the  constant  of  integration  has  been  ignored  for 
convenience)  . 


C C 

si  s, 

. 1 • 1 

s+c~!+m  s -0  • 

P e 

The  initial  conditions  on  s are  taken  to  be  gi.V'.n  by 
s(0)  - 0,  s (0)  = s 


(E-2) 


(E-3) 


where  sq  is  the  cage  velocity  at  the  onset  of  the  impact  between  the  cage 
and  the  ball.  The  form  of  the  solution  for  Equation  E-2  depends  on  the 
magnitude  of  the  parameter  4c2/m  C . If 

P 6 


D 

P 


4C 


MeC 


> 1 


(E-4) 


then,  the  solution  of  Equation  E-2  satisfying  the  initial  conditions  is 


so  -C2t 

s = — e sin  Cjt  , 


(E-5) 
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(a)  Ball-Cage  Impact 


F ■ C ( s-6) 


(b)  Force  Diagram 


FIGURE  E-l.  FORCE  DIAGRAM  FOR  BALL-CAGE  IMPACT 


where 


and 


The  rebound  velocity  s will  be  the  value  of  s when  C^t  - " or 


On  the  other  hand  if 


D 

P 


M C 


e s. 


< 1 


> 


(E-6) 


E-7) 


then  the  solution  is 


s "c2t 

9 = — 7 e sinh  c't  , (E-8) 

C1  1 


where 


(E-9) 


oscillatorvtnif)n°tef  ^at  the  solu,:lon  8lven  by  Equation  E-5  denotes  an 

ls  uistabll  ST  °f  Cafe'  ™1S  1S  taken  t0  ^ that  the  "Ob1”" 

s nstable . The  solution  given  by  Equation  E-8  indicates  that  the  ball- 

ST-SSVS1  abS°rb  the  e££eCtS  °f  the  Ca8e-ba11  bapact  an?  hence, 

£ "ge  ah™1dthbeCrbi:."°t  this  latter  condition,’ 


52 


APPENDIX  F 


MOMENTUM  TRANSFER  TO  CAGE -CONE  AT  TMPAflT 


APPENDIX  F 


MOMENTUM  TRANSFER  TO  CAGE-CONE  AT  IMPACT 


Ball-Pocket  Impact  for  Ball-Guided  Cape 


The  geometry  being  analyzed  is  illustrated  in  Figure  F-l  and 
involves  a cone-cylinder  cage  pocket.  Ball  contact  can  occur  either  at.  the 
cyiindrica1  surface  region  or  at  the  conical  surface  region.  The  impact 
conditions  will  be  assumed  to  be  similar  for  both  cases  though  somewhat 
more  complex  transformations  are  required  for  the  cone  contact. 


Description  of  Cone  Region 

. . _ /‘least  two  parameters,  in  addition  to  the  pocket  radius,  are 
nee  ed  to  define  the  location  of  the  cone  contact  relative  to  the  ball. 

(1)  The  cone  angle,  <5. 

(2)  The  perpendicular  distance  between  the  centei  of  the 

ball  and  the  base  of  the  cone  surface,  y . 

’ Js 

The  second  of  the  above  locating  coordinates  is  related  to  the  cage  design 
and  the  location  of  the  cage  center  of  mass  relative  to  the  ball.  That  is 


- y 


cm 


(F-l) 


where  yd  is  the  distance  between  the  center  of  the  ball  and  the  base  of 
the  cone  for  a perfectly  centered  cage  and  y is  the  radial  location  of  the 
cage  center  of  mass.  It  can  be  shown  that  cm 


Ycm  - P cos  (g  - ac)  , (F-2) 

^!re  p’  J*  aad  “c  are  the  locating  coordinates  for  a point  on  the  cage 
(Figure  3).  For  the  case  where  y = 0,  the  ball  will  make  simultaneous 
contact  with  the  cone  and  the  cylinder  if 


yd  = 


r 11  - c os  d) 
B sin  0 


(F-3) 


where  rg  is  the  ball  radius.  The  value  of  y 
represent  a reasonable  design  for  the  cage  to 


given  by  Equation  F-3  should 
be  ball  guided. 
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POCKET  CONTACT 


CONE  SURFACE 
v (TYPICAL) 


CONE  CONTACT 


FIGURE  F-l.  GEOMETRY 


Criterion  for  Contact- 


expressed 


II*  France  between  bell  and  cylinder  (Figure  F-l)  can  be 


Ar  ■ 


r 

P 


- r_  - 


(F-4) 


where  6 
ball. 


„ represents  the  tangential  movement  of  the 
rhe  corresponding  clearance  between  ball  and 


cage  relative  to  the 
cone  can  be  written 


Ac 


COS0 


+ yo  tan0  - 6 


(F-5) 


The  relationship  between  cone  and  cylinder  clearance 


can  be  expressed 


Ac  = Ar  - 


COS0 


c os  a v 

J+ 


tan0 


(F-6) 


It  can  be  seen  that  for  the  case  where  4 - n f 
F-3),  Ac  - Ar.  In  the  computer  routine’  it  is  I (“Sing  Equatlon 

of  Ac  versus  Ar  at  the  contact  conditions  to  esJablish^h^i,^6^  ^ leVel 
occurs  at  the  cylinder  or  at  the  cone  or  at  wl^1  J whether  contact 

ne  cone,  or  at  both  surfaces  simultaneously. 

Impact  and^  Rebound  Velocity  Components 

to  the  components  for  ball-cylinder°contact  but  only  E*gure  F_1  are  anal°8° 
transformation  of  the  ball-cvl  in/ior-  i t °nly  at  a skewed  angle.  Th 

velocities  can  be  Ztten  7 velocity  equation  into  cone  contact 


Jc  “ \ COS0  + sin0  , and 


(F-7) 


“ui  sin0  + v^  cos0 


(F-8) 


vr:L"%?  ‘Z  “T**  and  “i  and  V1  a«  th.  impact 
Appendix  D (Equation  D-4  and  D-9)  P Ml^r6TCt  V6ly’  discussed  in 
and  tangential  velocityZmpoLZ  ce^t^Z'’  ““  ”di*‘ 


u0  = uc  cos0  ■ vc  Sin0  , and 
o o 


(F-9) 


Vo  ~ uc  sin0  + v cos0 
o Co 


(F-10) 
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Here,  the  subscript,  c refers  to  the  rebound  condition  at  the  ball  cone 
interface. 


Equation  F-7  through  F-10  are  incorporated  into  the  transfor- 
mation subroutine  to  account  for  any  type  of  cone  contact.  If,  for  example, 
a zero  angle  is  used,  u^  and  v revert  to  simply  u.  and  v and,  concomitantly 

u and  v are  the  same  as  u and  v . 11 

o o c c 

o o 
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